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Abstract. Discriminating hadronic molecular and multi-quark states is a long standing problem in 
hadronic physics. We propose here to utilize relativistic heavy ion collisions to resolve this problem, 
as exotic hadron yields are expected to be strongly affected by their structures. Using the coalescence 
model, we find that the exotic hadron yield relative to the statistical model result is typically an order 
of magnitude smaller for a compact multi-quark state, and larger by a factor of two or more for a 
loosely bound hadronic molecule. We further find that some of the newly proposed heavy exotic 
states could be produced and realistically measured at RHIC and LHC. 
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We are now in a new stage of hadron physics, where various exotic hadron candidates 
have been discovered starting from a penta-quark state ® + (1530) [1] and tetra-quark 
states, D 4 y(2317) [2] and X(3872) [3]. We cannot properly explain these states within 
the simple quark model, then multi-quark and/or hadronic molecule components would 
be expected. An important aspect of exotic hadron physics thus involves the discrim- 
ination between a compact multi-quark configuration and a loosely bound molecular 
configuration. We have recently found that hadron yields in relativistic heavy ion col- 
lisions could provide useful information to address this question [4]. The hadron yield 
in the coalescence model relative to the statistical model result is found to be smaller in 
a compact multi-quark configuration because of the suppressed coalescence probability 
from additional quarks, and larger in a loosely bound hadronic molecule state because 
of the large size in which constituent hadrons can coalesce. 

In the statistical model, the number of produced hadrons of a given type h per unit 
rapidity is given by [5] 

^at =v jh_ r i*i (1) 

h 2k 2 . h r -i e (E h -B h n B -s h ns)/T H±1 > KJ 
where g%, B^, Sh are, respectively, the degeneracy, fugacity, baryon number and 



strangeness of the hadron. Following the expanding fire-cylinder model [6] for central 
Au+Au collisions at ^/snn = 200 GeV at RHIC, the volume, temperature, and baryon 
and strangeness chemical potentials of the source at statistical hadron production are 
taken to be V H = 1908 fm 3 , T H = 175 MeV, \x B = 20 MeV and Ms = 10 MeV, respec- 
tively. The fugacity y c = 6.40 is introduced for c and c quarks, to reproduce the expected 
total charm quark number N c = 3. 

In the coalescence model [6, 7], the hadron yield of type h containing n constituents 
at mid-rapidity is obtained using harmonic oscillator (Gaussian) wave functions as, 
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where Nj (gj) is the number (degeneracy) of the j'-th constituent, jli (/,■) is the reduced 
mass (orbital angular momentum) for the z'-th Jacobi coordinate, and o, ■= 1/ ^/JijCO 
with co being the oscillator frequency. Taking the quark numbers N u = Nj = 245 and 
N s = 150 at hadronization (T = 175 MeV, V = 1000 fm 3 ) for RHIC [6], we find that the 
addition of a s(p)-wave quark leads to the coalescence factor of 0.360 (0.093). Hadrons 
with more constituents or finite orbital angular momentum [8] are hence suppressed. 
By fitting the reference normal hadron yields (such as p, A(1115) or A c (2286)) in 
the statistical model, we fix CO = 550,519(385) MeV for hadrons composed of light 
quarks, and light and strange (charm) quarks, respectively. Weakly bound hadronic 
molecules are assumed to be formed at the kinetic freezeout point (7> = 125 MeV, 
Vf = 11322 fm 3 ). For a two-body s-wave hadronic molecule state, co is determined 
from the radius (co = 3/(2/i^(r 2 ))) or the binding energy (co = 6 x B.E.). 

We show the list of hadrons considered here in Table 1, including the proposed 
states, T c l c [9, 10], Z+(4430) [11], KKN[12], DN and DNN [13], @ cs [14], H [15], 



KNN[16], QCl [17], <mdH++ [10]. In Fig. 1, we show the ratio R h = N^ 1 / between 



the yields in the coalescence and statistical models. The grey band (0.2 < Rh < 2) 
covers the normal hadron ratios with 2q and 3q (open triangles), including the crypto- 
exotic hadrons with usual 2q/3q configurations. The ratio is below the normal band 
(Rh < 0.2) for a compact multi-quark configuration (diamonds). In particular, for light 
multi-quark configurations such as the tetraquark configurations of /o(980) and «o(980), 
the ratios are an order of magnitude smaller than those of normal hadrons or molecular 
configurations. This is consistent with a naive expectation that the probability to combine 
n-quarks into a compact region is suppressed as n increases. This suppression also 
applies to 5q states in multi-quark hadrons (A(1405), ©+(1530) and KKN) and the 
8g state in KNN. The ratios obtained by assuming hadronic molecular configurations 
(circles) are found to lie mostly above the normal band (Rh > 2). Moreover, these ratios 
depend on the size of the hadronic molecule; loosely bound extended molecules with 
larger size would be formed more abundantly. For A(1405) as an example, we find 
Rh = 1 • 1 for a small size (co = 174 MeV) obtained from the relation between the binding 
energy and CO, while a coupled channel analysis [18] gives a larger (r 2 ) and hence a small 
CO (= 20.5 MeV), leading thus to a larger R h = 4.9. 

We conclude from the above discussions that the yield of a hadron in relativistic heavy 
ion collisions reflects its structure and thus can be used as a new method to discriminate 
the different pictures for the structures of multi-quark hadrons. Indeed, based on the 



TABLE 1. Quantum numbers, configurations, and oscillator frequencies in hadronic molecule config- 
urations for exotics discussed in this proceedings. The undetermined quantum numbers, un-established 
particles, newly predicted particles, and p-wave coalescence configurations are marked by *\ $\ ^\ and 
(p) respectively. For hadron molecules, (BmoI. is fixed by the binding energy of hadrons (B), the inter- 
hadron distances (R), or the same as that for the subsystem (T). m and COmoI. are given in the unit of MeV. 
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ratio /o(980)/p ~ 0.2 from the preliminary measurement of the STAR Collaboration 
at RHIC [19] and the statistical model prediction p° = 42, we find the number /o(980) ~ 
8. This is enhanced from the statistical model result (~ 5.6). Comparing with the 
coalescence model results (3.8, 0.73, 0. 10, 13 for qq, ss, qqss, KK), the measured yield is 
consistent with the picture that the /o(980) has substantial KK components, and a pure 
tetraquark configuration can be ruled out for its structure. Further experimental effort to 
reduce the error is therefore highly desirable in understanding the structure of /o(980) 
and to put an end to this highly controversial issue. Similarly, efforts to measure the 
yields of other hadrons and newly proposed exotic states listed in Table 1 will provide 
new insights to a long standing challenge in hadron physics. 
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FIGURE 1. (Color online) Multi-quark hadron production at RHIC in the coalescence model relative to 
the statistical model. The patterns also holds for LHC as freezeout conditions are similar to that of RHIC. 
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